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Gravitational waves from a binary with a single dynamically significant spin, notably including 
processing black hole-neutron star (BH-NS) binaries, let us constrain that binary’s properties: the 
two masses and the dominant black hole spin. Based on a straightforward fourier transform of 
h{t) enabled by the corotating frame, we show the Fisher matrix for processing binaries can be 
well-approximated by an extremely simple semianalytic approximation. This approximation can be 
easily understood as a weighted average of independent information channels, each associated with 
one gravitational wave harmonic. Generalizing previous studies of nonprecessing binaries to include 
critical symmetry-breaking precession effects required to understand plausible astrophysical sources, 
our ansatz can be applied to address how well gravitational wave measurements can address a wide 
range of astrophysical and fundamental questions. Our approach provides a simple method to assess 
what parameters gravitational wave detectors can measure, how well, and why. 

PACS numbers: 


Introduction. - Ground-based instruments like LIGO 
[T] and Virgo [5]) will soon identify and measure the 
properties miH] of the relatively well-understood grav¬ 
itational wave (GW) signal from the nearly adiabatic 
and quasicircular inspiral of the lowest-mass coalescing 
compact binaries (CBCs) [T5H5D] : binaries consisting of 
either black holes or neutron stars with total masses 
M = mi -k m2 < I6M0 and intrinsic spins 81,82 that 
satisfy the Kerr limit |8i|/mf < 1. These measure¬ 
ments’ accuracy determines the range of astrophysical 
and fundamental questions that can be addressed via 
gravitational waves, including but not limited to identi¬ 
fying how coalescing compact binaries form [31I - I37j ; how 
the universe expands |38j ; how high-density nuclear mat¬ 
ter behaves and responds [35H1H]; and even how reliably 
general relativity describes the inspiral, coalescence, and 
gravitational radiation from each event [49]. In general, 
astrophysical formation channels [321 1501155] will popu¬ 
late generic spin orientations, not just high-symmetry, 
nonprecessing configurations with 81,82 parallel to L. 
For these most likely sources, spin-orbit and spin-spin 
couplings cause the misaligned angular momenta to pre- 
cess [50053] , breaking degeneracies present in the high- 
symmetry case and thus enabling higher-precision mea¬ 
surements [1 nni m ism]- While powerful analytic 
techniques were developed to estimate the measurement 
accuracy for nonprecessing binaries , then broadly 

applied [31037] [^ [30048] [49] . for processing binaries a 
comparable theoretical tool has remained unavailable. 
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Instead, the measurement accuracy has been evaluated 
on a case-by-case basis numerically, usually by Bayesian 
methods that systematically compare the data with all 
possible candidate signals [SIH]. 

The main result of this work is a generalization of the 
classic analytic approach used to approximate measure¬ 
ment accuracy [55] to the case of precessing binaries with 
a single significant spin undergoing an extended adia¬ 
batic, quasicircular inspiral. We restrict to a single signif¬ 
icant spin both for convenience - this limit is well-studied 
[5i[M]iin] - and without significant loss of generality - 
the smaller body’s spin often has little dynamically signif¬ 
icant impact on the angular momentum budget, orbit, or 
processional dynamics (I82I < <C |8i,|L|), allowing 
a single-spin model to adequately reproduce the dynam¬ 
ics and posterior m- To highlight the broad utility of 
our approach, we defer concrete but arbitrary implemen¬ 
tation details until we evaluate numerical results. Our 
result is important because it provides the first powerful 
analytic tool to assess what can be measured using grav¬ 
itational waves and why, includes the critical symmetry¬ 
breaking effects of spin precession. 

Inference from GW. - Bayes theorem provides an un¬ 
ambiguous expression for the posterior CBC parameter 
distribution given instrumental data {d}; see, e.g, PBOO 
m for a review. The signal and network response to a 
quasicircular GBC inspiral is characterized by eight in¬ 
trinsic parameters A [= (771.1,7712,81,82)] that uniquely 
specify the binary’s dynamics and seven extrinsic param¬ 
eters 9 [four spacetime coordinates and three Euler an¬ 
gles] that specify where, when, and with what orientation 
the coalescence occurred. Each detector responds to an 
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imposed strain as dk = rik + ReFk{'N)*h{t + N • x^jA, 0), 
where are the detector positions, Fk are antenna re¬ 
sponse functions, N is the line of sight to the source, 
(a member of 9), h{t\X, 9) is the strain derived from far- 
held solution to Einstein’s equations for the CBC; and 
rifc is some random realization of detector noise. The dis¬ 
tribution of stationary gaussian noise nk(t) in the A:th 
detector is completely characterized by its covariance 
or power spectrum (hfc(/)*hfc(/')) = ^Sk{\f\)S{f - /'). 
Let us dehne inner products on arbitrary com¬ 

plex functions a,b as (a|6)j, = 2 df s' 
log likelihood ratio C{9,9) favoring one signal with pa¬ 
rameters \,9 versus no signals is 21n£ = {dk\dk)k ~ 
{dk — hk{9, A)|dfc — hk{9, A)). If the data is known to con¬ 
tain a signal with parameters Aq = (Ao,0o)j we will de¬ 
note the parameters by A = (A, 9) and the likelihood by 
£(A|Ao,{n}). The signal amplitude p is set by the ex¬ 
pected value of the log likelihood {p^ = 21n£(A|Ao,0)). 
Using 15-dimensional posterior distribution ppost(A, 0) = 
Cp{9)p{\)/ J dXd9p{9)p{X)C{X,9), the measurement ac¬ 
curacy in some parameter Ai follows from the 90% con- 
hdence interval derived from a one-dimensional marginal 
distribution Ppost(Ai) = / d9dX2 . ■. dXsPpost{X, 9). 

While straightforward but expensive numerical tech¬ 
niques exist to estimate the marginal posterior distri¬ 
bution and hence measurement accuracy for concrete 
sources Aq and noise realizations nk [SHE], an equally 
straightforward analytic approximation to the (average) 
log likelihood exists at high signal amplitude, the Fisher 
matrix [HSHSHl [72U74] Tab. The Fisher information ma¬ 
trix arises in a quadratic-order approximation to the log- 
likelihood [ln£(A|Ao) ~ — irob(A—A*)a(A—A*)b-|-const, 
with A* the location of the noise-realization-dependent 
maximum]; often depends weakly on the specific noise 
realization used, particularly at high amplitude; and can 
be evaluated by a simple expression involving inner prod¬ 
ucts of derivatives. For example, for a source directly 
overhead a network with equal sensitivity to both polar¬ 
izations m, the Fisher matrix is 


/ dh dh \ 


( 1 ) 


The Fisher matrix can always be evaluated numerically 
via direct differentiation [52HS1]) henceforth denoted 
Fisher-D. For nonprecessing binaries, however, Poisson 
and Will [SS] introduced a powerful analytic technique, 
denoted here as Fisher-SPA: express the signal using 
a single dominant harmonic, with a necessarily-simple 
form; evaluate the fourier transform h via a stationary 
phase approximation; thereby evaluate the derivatives 
dah analytically] and, by reorganizing the necessary in¬ 
tegrals analytically, reduce the evaluation of Eq. 0 to 
an analytic expression and a handful of tabulated inte¬ 
grals. Despite its limitations, this method remains the 
most powerful and widely-used theoretical tool to esti¬ 
mate what can be measured and why. In the remainder 
of this work, we will review and generalize Fisher-SPA to 


processing binaries. 

As a matrix in at least 11 dimensions (15 with two 
processing spins), the Fisher matrix is both difficult 
to interpret and highly prone to numerical instability. 
For nonprecessing binaries, several studies have demon¬ 
strated that the intrinsic and extrinsic posterior distri¬ 
butions largely seperate [Til ESI and that the intrin¬ 
sic distribution depends weakly if at all on the spe¬ 
cific network geometry. Hence, to quantitatively as¬ 
sess what can be measured and why for a nonprecess¬ 
ing binary, it suffices to consider a source directly over¬ 
head and optimally aligned with a fiducial detector net¬ 
work [Ml ESHZZj- In the high-amplitude limit where 
the likelihood is well-approximated by a gaussian, the 
event time and polarization can be marginal¬ 
ized out analytically. Using this approximation, the 
log likelihood is approximated using an “overlap” P 

P(A, A') = maxt^ (h(A|£, V’c)|^(A|0,0))/||/i(A)||||/i(A')|| 
where ||/i|| = ^/Jh\ll), via Eq. (18) of [75]; the Fisher 
matrix arises as a quadratic approximation to the over¬ 
lap. This approach and its relatives, denoted here by 
Fisher-0 and in the literature by overlap, mismatch, or 
ambiguity function methods, has been widely adopted 
when analyzing numerical simulations |78H82j and cir¬ 
cumvents the numerical challenges that plague brute- 
force 11 -dimensional calculations. 

Even for processing binaries, several studies have sug¬ 
gested that the four spacetime coordinates decouple from 
intrinsic parameters and the binary’s three Euler angles 
[59l [83] . To simplify subsequent analytic calculations, 
following prior work IZSIITSI we will therefore adopt the 
ansatz that the source can be assumed directly overhead 
a network with equal sensitivity to both polarizations. 
Stationary-phase approximation. - The outgoing GW 
signal h{t, h, X) is modeled using a stationary phase ap¬ 
proximation of the leading-order (corotating) quadrupole 
emission, assuming a single significant spin. Specifically, 
adopting the conventions of [69| and PBOO, we express 
the strain for a source with intrinsic parameters A as a 
spin-weighted spherical harmonic expansion h{t\X,9) = 
h+ - ih^ = e~'^'^'^-^{M/dL)J2im^im{t - tc\X)Y‘^-^^{h), 
relative to a cartesian frame z' = J,x',y' defined by 
the total angular momentum J, where di is the lumi¬ 
nosity distance to the source, M is the binary mass, n 
is the propagation direction away from the source, £ is 
the coalescence time, and ipj is the angle of J on the 
plane of the sky. Within the post-Newtonian approxi¬ 
mation, both the amplitude and phase of these functions 
him {t) and their (stationary-phase) fourier transforms are 
slowy-varying and analytically-tractable functions [30j . 
For a nonprecessing binary, the sum is dominated by a 
single pair of complex-conjugate terms £22 = £2 -2 = 
1 ^ 22 1enabling efficient calculation of Eq. 0 [M] 
for an optimally-aligned source directly overhead a single 
detctor (Fisher-0). 

A generic quasicircular binary will orbit, precess, and 
inspiral on three well-seperated timescales l//orb; £recj 
and tinsp- For this reason, to a good approximation. 
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the gravitational radiation from a precessing, inspi- 
ralling binary |84L I85j can be approximated as if from 
an instantaneously nonprecessing binary: hi^(t\X) = 

where R{t) is a minimal ro¬ 
tation transforming z into L [SHI EOl [SB], where 
are available in the literature [lilSlISTj in terms of 
the spins, the orbital phase $ 0 ^ 6 , and a post-Newtonian 
expansion parameter v = {Md^orb/. The quanti¬ 
ties appearing in these expressions (L, S^, are 

determined by post-Newtonian approximations that pre¬ 
scribe the evolution of both spins [dtSi = fti x Si and 
the orbital phase ], [30l|88l|89|, where M is 

the rate at which the black holes’ mass changes jsniinii, 
henceforth neglected; R{v) is the rate at which energy is 
radiated to infinity; and E(v) is the energy of an instanta¬ 
neously quasicircular orbit; all of which are provided in 
the literature. At leading amplitude order, corotating- 
frame strain satisfies {t\X) = \him{t)\exTp{—i^orb), 
where \him\ is a slowly-varying function of u; substitut¬ 
ing this form into the general expression implies 

hUt\X) = (2) 

ffl 

In this expression, we have expanded the rotation R{t) 
using Euler angles, set by the orbital angular momen¬ 
tum direction L expanded relative to the (assumed 
fixed) total angular momentum direction J as L = 


sin (3jL cos a jlx' -I- sin/3ji sinaj^y' -|- cos/3ji,J; the re¬ 
maining Euler angle 7 = — J cos/3jLdajL- In this 
work, we restrict to the leading-order gravitational-wave 
quadrupole h^22 = — 81 / 71 / 577 ^^ exp(=F*$orb), so this 
sum has only two terms. 

For a binary with a single dynamically significant 
spin, the spin-precession equations imply that the or¬ 
bital angular momentum precesses simply around the 
total angular momentum: /3 changes slowly, on the in¬ 
spiral timescale, while a and 7 evolve on the precessional 
timescale [SBIIB^ . 

Because of seperation of timescales, because of the sim¬ 
ple form of Eq. ([^, and particularly because the phase 
terms m(<i> -I- 7 ) -b vna are monotonic and well-behaved, 
the stationary-phase approximation to the fourier trans¬ 
form = f dthimCxpicut can be carried out term 

by term [69) . For each term, the stationary-phase con¬ 
dition defines an m, m-dependent time-frequency trajec¬ 
tory Tmm(w) set by solving 


w = m(^orh 


d cos /3jl) + rnd 


(3) 


or, 

M 


equivalently, 


V 

-, 1/3 


m 


(M<i>;,ji/3 

Using this time-frequency 


relationship to evaluate '^rn,,fh = wt — m($orb + C) ~ 
and the slowly-varing coefficients in each term in Eq. 
([^, the stationary-phase approximation is 


hlm{uj) = 







•\/d™(*orb+C")+™“")/27r 


muj > 0 
fhui < 0 


(4) 


This expression for the SPA had been previously de¬ 
rived in the restricted PN approximation by precisely 
this method |69) . for simplicity neglecting the distinct 
time-frequency trajectories implied by Tm,Th¬ 
inner products via a time-frequency ansatz. - Most of 
the terms in Eq. ([^ are mutually orthogonal. For ex¬ 
ample, the modes with m > 0 and m < 0 have different 
helicity and have almost no overlap, indepdendent of the 
precession state [73 IMS]. Additionally, for binaries 
with more than a few precession cycles in band (see, e.g., 
[96j for suitable conditions on |Si|, mi, m 2 ), each term in 
Eq. <111 is associated with a unique time-frequency tra¬ 
jectory and hence is orthogonal to all others. Using this 
ansatz, the inner product {h{A)\h{A')) for A ~ A' can be 
approximated using a sum over 10 terms: 5 for each of 


the I = 2 modes m = — 2 , — 1 ,... 2 and, for each mode, 
one term for each helicity. 

By contrast, to evaluate the overlap of terms that are 
not orthogonal, the specific time-frequency trajectory 
has relatively little impact. We therefore approximate 
Tmrh — Torn henceforth. 

Fisher matrix. - We now use the time-frequency ansatz 
and the restricted PN approximation to evaluate a Fisher 
matrix for a source directly overhead an idealized net¬ 
work of two interferometers oriented to have equal sensi¬ 
tivity to both polarizations [Eq. Qj. Because of the 
time-frequency ansatz, the overall Fisher matrix is a 
weighted sum of 10 individual Fisher matricies, associ¬ 
ated with each harmonic: 
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2 

= E E P^^rns^Tb 

m—2 s—± 1 

2 = \j2 (o\l2 f°° d.f (~r M 

P2ms — l^2m J ^ )^m,2s\l^)\ J Sfi{f) SdP' V^-^c/j 

^(ms) /o°° ;sf^(7r7Wc/)"^/^i9a(«'2-2C-msa)ab(«'2-2C-msa) 


(5) 

( 6 ) 
(7) 



j A 

FIG. 1: Comparison: Comparison between marginalized 
Fisher matrix in 7V4c,?7 for a precessing binary evaluated us¬ 
ing the ansatz in this paper (red solid curve: Fisher-SPA); via 
full 15-dimensional MCMC (thin curves); and via a quadratic 
approximation to the 7-dimensional overlap (Fisher-O). The 
latter two results, previously presented O’Shaughnessy et al. 
[59], used with permission. 

In this expression, 1112 = uit — 2^orb is the stationary- 
phase approximation derived via ui = 2^orb- The weights 
Pms are associated with the relative contribution of each 
model m and sign s to the detected amplitude, along this 
line of sight. The 10 individual Fisher matrices Fq;, re¬ 
flect the Fisher matrix implied by a single harmonic, with 
a modified phase versus time to reflect that harmonic’s 
precession-induced secular phase change; each one re¬ 
duces trivially to the well-known nonprecessing Fisher 
matrix in the absence of precession. 

Because measurements often cannot tightly constrain 
all parameters, whether computed directly (Fisher-O) or 
via our approximation, the Fisher matrix is often degen¬ 
erate, particularly in phase angles. Following prior work, 
when evaluating the Fisher matrix to produce our final 
numerical results, we adopt a very weak (regularizing) 
prior to break degeneracy: F final =T -\- K where K re¬ 
flects a multivariate gaussian distribution with standard 
deviation 27r in phase angles, IMq in total mass, 1 second 
in time, and 1/4 in mass ratio. 

Implementation and results. - Figure [l] shows a compar¬ 


ison beteween our approximation; a Fisher-O approxima¬ 
tion; and a full MCMC posterior distribution. In our cal¬ 
culation, to minimize superfluous differences associated 
with uncontrolled post-Newtonian remainders, we eval¬ 
uated the phase functions 4>2(u), 7 (u),a(u) needed for 
our Fisher-SPA approximation numerically, using pre¬ 
cisely the same post-Newtonian evolution model adopted 
in the other calculations shown (SU] to evaluate t{v) 
and $ 2 ( 1 ),a(t),7(t); see that work for the specific post- 
Newtonian approximation used. Despite not including 
the prior p(A) and despite adopting highly simplifying 
assumptions, our expression shows good agreement with 
the multidimensional posterior and previous numerical 
estimates of the marginalized Fisher-O matrix. Though 
not shown here, similar results follow by using the explicit 
expressions for 'I' 2 (u), 7 (u),a(z;) available the literature 

m- 

Future directions. - While only approximating the re¬ 
sults of detailed Bayesian parameter estimation [BIlZl 
IU5] . Fisher-matrix calculations provide a powerful and 
analytically-tractable tool to assess what can be mea¬ 
sured and why. Extending the Fisher-SPA method to in¬ 
clude a single precessing spin will help rapidly interpret 
of real gravitational wave data, via improved methods to 
explore the model space and interpret the posterior; as¬ 
sess the impact of systematic errors from the waveform 
model; quantify the accuracy to which tidal effects and 
modifications of general relativity can be detected; and 
otherwise understand what can be measured and why. 

Further investigation is needed to generalize our ap¬ 
proach to account for a second significant spin, using 
recently-developed analytic solutions [HISTII^; and 
to perform a large-scale comparison between our calcu¬ 
lations and detailed Bayesian parameter estimates. Fi¬ 
nally, to facilitate the immediate use of our approach and 
enhance its similarity to prior work, we have adopted ex¬ 
tremely simple assumptions (e.g., the neglect of r^,™; 
the neglect of all but 10 overlaps; and the restricted PN 
expansion). These assumptions can easily be relaxed if 
more detailed calculations are required, since the simple 
form of Eq. ([^ insures a theoretically-tractable analysis. 
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